EPJ manuscript No. 

(will be inserted by the editor) 



BRAG2007 Workshop SUMMARY 



S '■ The physical meaning of scattering matrix singularities 
o ' in coupled-channel formalisms 

(N : 

>■ '. S. Capstick\ A. Svarc^, L. Tiator^, J. Gegelia^, M. M. Giannini"*, E. Santopinto**, C. Hanhart^, S. Scherer^ 
O , T. -S. H. Lee6'7, T. Sato^'^ and N. Suzuki^'*^ 

^ Department of Physics, Florida State University, Tallahassee, FL 32306-4350, USA 
Rudjer Boskovic Institute, Bijenicka cesta 54, P.O. Box 180, 10002 Zagreb, Croatia 
' ^ Institut fiir Kernphysik, Johannes Gutenberg-Universitat Mainz, Johann-Joachim-Becher-Weg 45, D-55099 Mainz 
' * Dipartimento di Fisica, Universita degli Studi di Genova & I. N.F.N. Sezione di Genova 
(-H , ^ Institut fiir Kernphysik, Forschungszentrum Jiihch, 52425 Jiilich, Germany 
O ^ Physics Division, Argonne National Laboratory, Argonne, IL 60439 

I ^ Excited Baryon Analysis Center, Thomas Jefferson National Accelerator Facility, Newport News, Va. 22901 
Department of Physics, Osaka University, Toyonaka, Osaka 560-0043, Japan 

^ . 

Received: date / Revised version: date 

Abstract. The physical meaning of bare and dressed scattering matrix singularities has been investigated. 
Special attention has been attributed to the role of well known invariance of scattering matrix with respect 
to the field transformation of the effective Lagrangian. Examples of evaluating bare and dressed quantities 
in various models are given. 

PACS. ll.80.Gw, 13.40.Gp, 13.60.Le, 13.85.Fb, 13.75.Cs, 14.20.Gk, 14.40 Aq, 25.20.Lj, 25.40.2h 

This paper is a collection of different, sometimes con- to be made. One either has to dress quark-model reso- 
flicting standpoints presented at the BRAG2007 pre-meeting nant states spectrum and compare the outcome to the 
of the NSTAR2007 Workshop. It is organized in the fol- experimental scattering matrix poles, or to try to take 
lowing way: into account all self-energy contributions which are im- 

; I ■ The Introduction is written by A. Svarc, Sect. 2.1 by S. plicitly included in the measured scattering matrix pole 
^ Capstick, 2.2 by C. Hanhart, 3.1 by S. Scherer, 3.2 by J. parameters, make a model independent undressing proce- 
Gegelia, 4.1 by M. Giannini and E. Santopinto, 4.2 by T. - dure and compare the outcome to the impulse approxima- 
S. H. Lee, T. Sato and N. Suzuki and the Conclusion by tion quark-model calculations. The first options seems to 
A. Svarc, S. Capstick and L. Tiator. be feasible but complicated, but the latter one seems to be 

impossible due to very general field-theory considerations. 

^ , . , . . We report on investigating both options. 

1 Introduction b b t- 

An attempt how to un-quench the constituent quark 

Establishing a well defined point of comparison between model of ref . [T] , together with describing all accompany- 

experimental results and theoretical predictions has for ing complications, is presented. The procedure seems to 

decades been one of the main issues in hadron spectroscopy, be cumbersome, but straightforward. 

and the present status is still not satisfactory. Experi- „i i , • i • ; i ■ . n i 

, . ,. 1 m-iirA \ J Ti J 1 ■ i he second option, undressmg the experimentally ob- 

ments, via partial wave PWA and amplitude analysis , . , ,, • , • • , •,• i ; 

. T 1 1 • r ■ 1 ■ tamed scattering matrix singularities, however, seems to 

(AA), can give reliable information on scattering matrix , . , i i i i i i ; i 

, 1 1 1 1 i- 11 • be inherently model dependent due to very general ar- 

smgularities, while quark model calculations usually give , .. ,. ^ ^ , . . , ^ i i ; i a ■ ^ 

. r .. i i i i • ii £ i J guments originating from the local field theory. A simple 

information on resonant states spectrum m the first order , , .„ , ; . t, . , . . , i 

1 ■ /I / 1 1 j_ \ model, illustrating this claim is presented, 

impulse approximation (bare/ quenched mass spectrum). " 

And these two quantities are by no means the same. Up In spite of looking entirely dissimilar, problem of model 
to now, in the absence of a better recipe, these quantities independent undressing of full scattering matrix singular- 
have usually been directly compared, but the awareness ities seems to be strongly correlated to the recent con- 
has ripen that the clear distinction between the two has troversy whether the off-shell effects are measurable or 
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not. It is therefore essential to extend the discussion to 
(un)measurability of off-shell effects as well. 

For decades the off-shell properties of two-body ampli- 
tudes seemed to be a legitimate measurable quantity, and 
numerous attempts to get hold of it in nuclcon-nuclcon 
bremsstrahlung and real and virtual Compton scattering 
on the nucleon have been made. However, in early 2000-es 
it became apparent that strong field-theoretical arguments 
do not speak in favor of this claim It seems that 

it is very likely that the well known invariance of scat- 
tering matrix with respect to the field transformation of 
the effective Lagrangian [S] makes it possible to transform 
the off-shell effects into the contact terms for diagrams of 
the same power counting level. This effectively makes the 
off-shell effects an unmeasurable quantity. 

When applied to the effective two body meson-nucleon 
amplitudes, this statement implies that the ability of cou- 
pled-channel formalisms to separate the self-energy term 
and evaluate the bare scattering matrix poles (singulari- 
ties in which the meson-exchange effects are fully taken 
into account) is a model dependent procedure. Namely, 
any method for evaluating self-energy contributions un- 
avoidably demands a definite assumption on the analytic 
form of the off-shell interaction terms, hence introduces 
model dependent and consequently unmeasurable hadronic 
shifts. 

The invariance of different parameterizations of the 
scattering matrix singularities with respect to field - re- 
definitions is also the object of our study. Scattering ma- 
trix poles are nowadays quantified in two dominant ways: 
either as Breit-Wigner parameters, i.e. parameters of a 
Breit-Wigner function which is used to locally represent 
the experimentally obtainable T-matrix, or as scattering 
matrix poles (either T or K). In spite of the fact that it 
is since Hoehler's analysis [61I3 quite commonly accepted 
that Breit-Wigner parameters are necessarily model de- 
pendent quantities, they are still widely used to quantify 
the scattering matrix poles. Only recently the scattering 
matrix poles are being shown in addition. We demonstrate 
that within the framework of effective field theory, scat- 
tering matrix poles are, contrary to Breit-Wigner param- 
eters, unique with respect to arbitrary field-redefinitions. 

Bare and dressed scattering matrix quantities have for 
more then a decade been calculated and presented within 
a framework of various coupled-channel models [8lf9] , and 
definite correlation between scattering matrix singulari- 
ties and quark-model quantities has been in general es- 
tablished [10]. However, the most direct connection be- 
tween full scattering matrix singularities and hadron mod- 
els with confinement forces has been offered in [TT l lT2 t 
113) for the various versions of dynamical coupled-channel 
model. In these models the bare N* states are understood 
as the excited states of the nucleon if its coupling with 
the reaction channels is turned off, so the authors natu- 
rally speculate that the bare N* states of these models 
correspond to the predictions from a hadron model with 
confinement force, such as the well-developed constituent 
quark model with gluon-exchange interactions. The role 



of hadronic shift model dependence, however, is not ex- 
plicitly discussed. 

A simple conclusion emerges: dressed scattering ma- 
trix singularities are the best, model independent meeting 
point between quark model predictions and experiments, 
and bare quantities in coupled-channel models remain to 
be legitimate quantities to be extracted only within a 
framework of a well defined model. To understand and 
interpret them correctly, one has to keep track of the ex- 
istence of the hadronic mass shifts produced by off-shell- 
ambiguities, and take them fully into account. 



2 Dressing and undressing scattering matrix 
singularities 

2.1 Un-quenching the quark model 

The usual prescription for calculation of the masses of 
baryons is to ignore the effects of decay-channel couplings, 
which is the assumption that the states are infinitely long- 
lived. Given that baryon widths are comparable to the 
mass splittings between similar states caused by short- 
range interactions between the quarks, the effects on baryon 
masses of continuum (baryon-meson) states, or equiva- 
lently qqq-qq components, clearly cannot be ignored. The 
problem is that there are many distinct intermediate states 
which can contribute substantially to the self energies of 
baryons through baryon-meson loops, because of the pres- 
ence of many thresholds in the resonance region. Calcu- 
lations of the effect of two-meson intermediate states in 
mesons have been carried out, especially for the interesting 
problem of the cu-p mass difference [l4l[T5] which illustrate 
the complications which arise in the case of baryons. In or- 
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Fig. 1. Calculation of baryon self energies in the quark model 



der to calculate the self energy of a baryon -8(0) due to a 
particular baryon-meson intermediate state i3'(— k)C(k), 
as in Fig. [l] we require a calculation of the dependence 
of the vertex M.BB'c{k) on the magnitude k of the loop 
momentum k. This in turn requires a model of the spec- 
trum (including states not seen in experiment), which pro- 
vides wave functions for the baryons, and a model of the 
-6(0) -B'(— k)C(k) decay vertices. A popular choice for 
the former is some form of constituent quark model, and 
for the latter is a pair-creation model such as the ^Pq 
model illustrated in Fig. [21 where baryons decay by the 
creation of a quark-antiquark pair with the quantum num- 
bers of the vacuum. 
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Fig. 2. Pair creation model of baryon decays 



In order to self-consistently calculate the masses of 
baryons in the presence of baryon- meson intermediate states, 
one possible approach [16] is as follows. The masses and 
decays are calculated using a three-quark Hamiltonian 
Hqqq and a pair-creation Hamiltonian i?pc, that depend 
on strong coupling, quark mass, and string tension pa- 
rameters a1, TO?, and 6", etc., and a pair-creation coupling 
strength 7. These parameters are usually determined by a 
fit to the (dressed) spectrum Eb and decay partial widths 
in the absence of Fock-space components higher than qqq, 
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The correction due to the loop B ^ B'C \s 
Eb - M^(a°, mO, 6°, ...) + Sb'^Eb-Eb'^uI, to^, 6"), 



where 



^B'C 



|(i?'(-k)C7(k)|i/pe|i?(0))r 



Eb ~ y/El^+P - y/El+¥ + 



The imaginary part of the loop integral is Fb-^b'C l"^- A 
sum is to be performed over baryon-meson intermediate 
states B'C , and the parameters a°, to°, 6°, and 7,... are 
to be adjusted for self-consistent solution with Eb equal 
to the dressed baryon mass. In principle one should solve 
similar equations for the meson masses Ec- 

This procedure is equivalent to second order perturba- 
tion theory in the decay Hamiltonian -ffpc, and allows cal- 
culation of the (momentum-space) continuum B' C com- 
ponent of the dressed baryon states, and also the mixing 
B — >• B'C — > B" between different baryon states caused 
by the continuum intermediate states. 

Calculation of this kind have been applied to N , Z\, A, 
S and S* ground and (singly) orbitally excited states us- 
ing intermediate states made up of ground state baryons, 
with the pseudoscalar mesons tt, K, rj, 77' in Ref. [17] 
and these pseudoscalar mesons plus the vector mesons p, 
oj and K* in Refs. [18,19,20J. Because there are many 
baryon-meson thresholds nearby in energy, for example 
the Np and Ap thresholds are close to those of iV(1535)7r 
or A{1A05)K , one should not restrict the intermediate me- 
son states to TT, or even all pseudoscalars, or the interme- 
diate baryon states to N and A, or even all octet and 
decuplet ground states. 

Zenczykowski fW showed that if one assumes exact 
SU(3)/(8)SU(2)spin symmetry, that only ground state ba- 
ryons and mesons exist, that all octet and decuplet baryons 
have the same mass and the same wave function, and 
that all pseudoscalar and vector ground-state mesons have 
the same mass Mq and the same wave function, that 



all self energy loop integrals are the same, apart from 
SU(6)spin-flavor factors at the vertices. Under these con- 
ditions we expect the sum of self energy contributions to 
the nucleon and Z\(1232) masses to be identical. Interest- 
ingly, the sum of the squares of the SU(6)spin-fiavor factors 
is the same only if we include all baryon-meson combina- 
tions (non-strange, strange, or both) consistent with the 
conserved quantum numbers, including both pseudoscalar 
and vector mesons. This is true of the self energies of any 
ground state baryon, and is also true if the ^Pq model 
is used to calculate the vertex factors, as it reduces to 
SV{6)w in this limit. 

Away from the SU(3)/ limit, Tornqvist and Zenczy- 
kowski [21] were able to show that with the introduc- 
tion of a simple pattern of SU(3) / breaking present in 
the ground-state baryon and meson mass spectra, that 
the usual SU(6) relations for baryon masses are present 
in the dressed baryon masses calculated to first order in 
the symmetry breaking parameters. This suggests that we 
can interpret SU(6) symmetry breaking effects as partly 
due to spin and flavor-dependent interactions between the 
quarks, and partly due to loop effects. 

It is clear from this and other calculations that the 
effects of these self energies on the spectrum are substan- 
tial. Zenczykowski I8J finds many mass splittings close to 
those of the dressed pole parameters from analyses, with- 
out spin and fiavor-dependent interactions between the 
quarks. Other calculations show splittings in the dressed 
P-wave (lowest orbitally) excited baryons which resemble 
spin-orbit effects [I7 i ri9 t ,20j. These could cancel against 
those expected from other sources and provide a solu- 
tion to the spin-orbit problem in certain quark models 
of baryon masses. 

These calculations lack a self-consistent treatment of 
external and intermediate baryon states, and so it is not 
clear that the sum over intermediate baryon-meson states 
has converged. Geiger and Isgur [TJ demonstrated that 
this sum does converge, albeit slowly, using a non-relativistic 
quark model for baryon masses and wave functions and 
a ^Pq model for their decays. Using a covariant model 
based on the Schwinger-Dyson Bethe-Salpeter approach 
was shown to lead to faster convergence in Ref. [TS]. A 
study of the dressed masses of N and A ground and P- 
wave excited baryons [22j which involves intermediate pseu- 
doscalar and vector ground-state mesons and many in- 
termediate baryons (ground states, and N, A, A, S, S* 
excited states up to the second band of negative-parity 
states at roughly 2100 MeV), representing hundreds of in- 
termediate states, is underway. Vertex form factors are 
calculated analytically using mixed relativized-model |23| 
wave functions and the ^Pq model [24] . 

This study shows that the usual ^Pq niodel gives ver- 
tices which are too hard, giving large contributions from 
high loop momenta. They can be softened by adopting a 
pair-creation form factor which decreases as the relative 
momentum of the created quark and antiquark increases. 
This calculation is currently being reworked to allow self- 
consistent renormalization of the quark model parameters. 
As an example, in Ref. [22l the strong coupling parameter 
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was reduced in order to take into account the additional 
A — N splitting in the sums over baryon-meson loops con- 
tributing to the self energies of both of these states. Simi- 
larly, in their calculation of these effects in mesons, Geiger 
and Isgur [Hj showed that the formation of an interme- 
diate meson pair was equivalent to string breaking, which 
has the effect of renormalizing the meson string tension. 
Barnes and Swanson 25J have examined shifts in the char- 
monium spectrum due to D, D* , Dg and D* meson pairs. 

In conclusion, the next Fock-space component is likely 
more important than differences among qqq models. Cal- 
culating its effects requires the use of a full set of SU(6)- 
related intermediate states, spatially excited intermediate 
baryons, and a careful treatment of mixing effects. Renor- 
malization of the parameters in the quark model parame- 
ters such as as, the quark masses, the string tension, and 
the pair-creation strength needs to be carried out system- 
atically. This requires examining the mass shifts of more 
than just iV, A and their negative-parity excitations. De- 
cay vertices need additional suppression when the dressed 
masses of external states are well above the threshold for 
an intermediate state, which is the case in relativistic mod- 
els. 



2.2 Undressing the dressed scattering matrix 
singularities 

To get a better understanding of the relation of bare quan- 
tities to dressed quantities it is sufficient to study a system 
of two nucleon like states [N and R) coupled to a scalar 
field [a) [26J. A possible Lagrangian reads 

Ci=N {i^ - Mn) N + R{i^^ Al'^) R 
+ ^{d^'adf,a^m^)+ga{RN + NR)... . (f) 

Here the superscript indicates that masses are bare quan- 
tities that undergo dressing beyond tree levefl The result- 
ing vertex is shown as diagram (a) in Fig. [31 The dots indi- 
cate possible more complex terms, like contact terms of the 
type a^RR (see Fig. [3l[c)). However, in phenomenological 
studies those are rarely included. From this Lagrangian we 
may now calculate observables like scattering amplitudes. 
To keep things simple we focus only on the self-energy of 
the R field. The corresponding diagram is shown in Fig. 
|4ja). The real part of this diagram provides the so called 
hadronic shift — the difference between the bare mass and 
the physical mass — and the imaginary part the width. 

A theorem based on very general assumptions in field 
theory states that if two fields (p and x related non- 
linearly ((j) — xF{x) with F{0) — I) then the same observ- 
ables arise if one calculates with using C{4>) or with x 
using C{xF{x)) [27|. Thus instead of the fields in Eq. ^ 
we may switch to a modified nucleon field defined through 



N 



N' = N + aaR , 



^ In principle also the coupling g and the mass of a and 
A'" are bare quantities, however, to ease notation we drop the 
corresponding superscript, for in what follows we focus solely 
on the self energy of the R field. 



N R N R R 

(a) (b) 



(c) 



Fig. 3. Vertices from the interaction Lagrangians of Eqs. ([T]) 
and ©. 



where a is an arbitrary, real parameter. Expressed in terms 
of N' the interaction part of the original Lagrangian now 
reads 

C{, = ga {RN' + N'R) ~ aa {R {i^ - Mn) N' + h.c.) 
+ aR[a^ {i^- Mm) -2ga]Ra + ... . (2) 

In addition to the vertex of the previous Lagrangian now 
two new structures appear: a momentum dependent RNa 
vertex (depicted in Fig. [Htb)) and a a^RR vertex (de- 
picted in Fig. [Sljc)). The resulting contributions to the 
i?-selfenergy from this Lagrangian are shown in Fig.jja)- 
(d). The field theoretic theorem quoted above gives that 
the total self energy contribution from the modified La- 
grangian is identical to that of the original Lagrangian 
with the same parameters. Especially, the hadronic shift 
remains the same and it seems that indeed it is a well 
defined quantity. However, the problem is that we do not 
know the true hadronic Lagrangian. Thus, starting from 
Eq. (Ill) is as justified as starting from the following inter- 
action Lagrangian 



Ci ^ ga {RN + NR) - aa (i? {i^ - Mn) N + h.c.) 

+ ... . (3) 



Obviously the only difference to the previous equation is 
that the a^RR vertex was abandoned. On the one loop 
level thus the only difference compared to the previous 
expression for the i?-selfenergy is that tadpole diagrams 
were removed. Since this class of diagrams does not lead to 
non-analyticities, their effect can always be absorbed into 
the bare mass and the wave function renormalization of 
the R field. Therefore, with properly adjusted parameters, 
the self-energy is the same to one-loop between the theory 
that follows from Ci and that from C2. 

Is there any way in practise to decide, which one of 
the two Lagrangians is to be preferred? The answer to 
this question is no for the following reasons: although the 
a^RR contact term can contribute to the R self energy 
at three loop order, this is of no practical significance, 
since not only has any effective Lagrangian a too limited 
range of applicability and accuracy to allow for the extrac- 
tion of such effects but also a complete treatment should 
include anyway direct aR —>■ aR transitions in both La- 
grangians in addition to the terms given explicitly above. 
The latter argument also applies to information deduced 
from aR —>■ aR cross sections. Therefore there is in prac- 
tice no way to decide which one of the two interaction 
Lagrangians — Eq. ^ or Eq. ^ — is to be preferred. As 
outlined above, however, quantities like the self-energies 
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(a) (b) (c) (d) 

Fig. 4. Selfenergies for the R field to one loop order from toy-model I (Eq. ([T}) and toy-model II (Eq. |3]). 



of the resonance R are different in the two approaches 
and consequently the bare masses as extracted from fits 
to experiment are different. We therefore conclude that 
bare masses (or in general bare quantities) do not have 
any physical significance. 

The question studied here is very closely linked to the 
question of measurability of off-shell effects. The argu- 
ment just presented can also be used as yet another il- 
lustration that off-shell effects are not observable. This 
already follows from a comparison of Eq. ^ and Eq. ^ . 
As argued above both lead to identical observables. Es- 
pecially, the on-shell RNa vertex that can be related to 
the decay width from R aN is the same for both mod- 
els. However, in our example for off-shell nucleons the 
vertex can be anything. In general, within a consistent 
field theory off-shell effects either can be absorbed into 
counter-terms or have to cancel exactly. The same issue 
is discussed for bremsstrahlung in Ref. [4j. Another illus- 
trative example of the cancellation of off-shell effects is 
provided in Ref. [.28J for the reaction NN NNn. 

It should be stressed that the question in focus here is 
very different to that of the relation between two-nucleon 
and three-nucleon observables and the presence of three- 
body forces. The main difference is that in the few nu- 
cleon systems it is possible to construct three-body forces 
that are consistent with the two-nucleon interaction used, 
e.g. within effective field theory — for a recent review 
see Ref. . Changing the two-nucleon interactions leads 
also to controlled changes in the three-body forces in the 
sense sketched above. However, what would be needed for 
a model independent extraction of bare hadron masses 
would be a method to identify the hadronic interaction 
that is the one that matches to the particular quark model, 
thus a connection is needed between two systems with 
very different degrees of freedom. We argue that it follows 
from the reasoning above that this identification can not 
be made as a matter of principle. However, the inclusion 
of hadronic loops within the quark model, as sketched, 
e.g. in the presentation by Simon Capstick, is obviously 
justified. 

We therefore have to conclude that the only quantities 
relevant for spectroscopy that can be extracted from ex- 
periment are resonance poles and the corresponding residues 
However, this is still a lot for both quantities contain im- 
portant structure information like the amount of SU(3) vi- 
olation or even the very nature of the state [30] . A extrac- 
tion of poles and residues from the data needs coupled- 
channel codes of the type of Refs. [5T1I5^1[55] with the cor- 
rect analytical properties and consistent with unitarity. 



Only the then a controlled analytic continuation to the 
complex plain is possible. 



3 Field-theory considerations 

3.1 From ofF-shell to on-shell kinematics 

It is a natural and legitimate question to ask whether 
the off-shell behavior of particular interaction vertices is 
unique and whether it is possible to extract such behavior 
from empirical information similarly as one, say, extracts 
the electromagnetic form factors of the nucleon from elas- 
tic electron scattering. In this context one might think of 
the electromagnetic interaction of a bound (off-shell) nu- 
cleon or the investigation of the off-shell nucleon-nucleon 
amplitude entering the nucleon-nucleon-bremsstrahlung 
process. For the case of pions, Compton scattering [34] 
and pion-pion bremsstrahlung [2] were discussed using chi- 
ral perturbation theory (ChPT) at lowest order. It was 
shown that off-shell effects with respect to the effective 
pion fields depend on both the model used and the choice 
of representation for the fields. From that the conclusion 
was drawn that off-shell effects are not only model depen- 
dent but also representation dependent, making a unique 
extraction of off-shell effects impossible. The spin-1/2 case 
was discussed in Ref. [3 . 

A related situation occurs when one is interested in 
corrections to current-algebra results obtained from the 
partially conserved axial-vector current (PCAC) relation 



(4) 



where A'^-^ is the isovector axial- vector current and is 
a renormalized field operator creating and destroying pi- 
ons; and -FV denote the pion mass and decay constant, 
respectively. While predictions of current algebra and the 
PCAC relation involve the so-called soft-pion limit, 
limqQ_>o limg^o[' ' amplitudes for physical pions are to 
be taken at — M^. Can the connection between soft- 
pion kinematics and on-shell kinematics be uniquely de- 
termined? The answer is yes, if the problem is entirely 
formulated in terms of the relevant QCD Green functions. 

We will illustrate these issues in the framework of ChPT 
[55] which establishes a systematic connection with the 
underlying field theory, namely, QCD. Let us first discuss 
off-shell effects with respect to the effective fields. To that 
end, we consider tttt scattering at lowest order in ChPT 
(see Sect. 4.6.2 of Ref. [36] for more details): 



Tr [df.Uid^Uy 



F^A'P 

££i!fzLTr(C/ - 



6 



S. Capstick et el: The physical meaning of scattering matrix singularities 



where — 2Brh. B is related to the quark condensate 
{qq)a in the chiral limit and rh is the average of the u- and 
d-quark masses 1351 : U is an SU(2) matrix containing the 
pion fields. We will use two alternative parameterizations 
of [/: 

1 



Uix) = — [aix) + iT ■ n{x)] , a{x) = ^F^^tt^x), 



U{x) 



cxp 



F^ 



) sm 



F., 



The a and exponential parameterizations are related by a 
field transformation (change of variables) 



TT 



1 



1 02 

6f! 



The relevant tttt interaction Lagrangians read 



1 A/2 



24:Fr^ ' 



Observe that the two interaction Lagrangians depend dif- 
ferently on the respective pion fields. For Cartesian isospin 
indices a, 6, c, d the Feynman rules for the scattering pro- 
cess 7r°(pa) + 7r*(pfc) — > 7r'^(pc) + T^'^iPd) read, respectively, 



F2 



F2 



F2 



■^2^ = -Mf 



3F2 



where we introduced Ak = pi ~ and the usual Man- 
delstam variables s = (pa + pbY , t — {pa —pcY, and 
u = {pa- PdY satisfying s + t + u = pi + pi + pi + p^. 
If the initial and final pions are all on the mass shell, i.e., 
Ak — 0, the scattering amplitudes are the same, in agree- 
ment with the equivalence theorem of field theory [37. 
On the other hand, if one of the momenta of the external 
lines is off mass shell, the amplitudes A^l'^ and AA'^'^ dif- 
fer. In other words, a "direct" calculation of tttt scattering 
in terms of the effective fields gives a unique result inde- 
pendent of the parameterization of U only for the on-shell 
matrix element. 

According to the standard argument in nuclcon-nuclcon 
bremsstrahlung one would now try to discriminate be- 
tween different on-shell equivalent tttt amplitudes through 
an investigation of the reaction TT°'{pa)+T^^{Pb) t^''{Pc) + 
'^'^{Pd) +l(k). This claim was critically examined in Refs. 
[2113]. To that end the electromagnetic field is included 
through the covariant derivative -D^C/ = 9^t7 +ieA^[Q, U] 
where Q = diag(2/3, —1/3) is the quark charge matrix. In 
the fj parameterization, the total bremsstrahlung ampli- 
tude is given by the sum of only such diagrams, where 
the photon is radiated off the initial and final charged pi- 
ons, respectively. One may then ask how the different off- 
shell behavior of the tttt amplitude of M.'^'^ enters into the 



calculation of the bremsstrahlung amplitude. Observe, in 
this context, that the exponential parameterization gener- 
ates electromagnetic interactions involving 2n pion fields, 
where n is a positive integer. In the exponential parame- 
terization an additional 4(^7 interaction term relevant to 
the bremsstrahlung process is generated. Hence the total 
tree-level amplitude now contains an additional four-pion- 
one-photon contact diagram. Combining the contribution 
due to the off-shell behavior in the tttt amplitude AA"^ 
with the contact-term contribution, we found a precise 
cancelation of off-shell effects and contact interaction such 
that the final results are the same for both parameteriza- 
tions. This is once again a manifestation of the equiva- 
lence theorem |37j . What is even more important in the 
present context is the observation that the two mecha- 
nisms, i.e. contact term vs. off-shell effects, are indistin- 
guishable since they lead to the same measurable ampli- 
tude. 

Now, what about the off-shell behavior of QCD Green 
functions? The method developed by Gasser and Leutwyler 
|35j deals with Green functions of color-neutral, Hermi- 
tian quadratic forms involving the light-quark fields q — 
{u, d)^ of QCD and their interrelations as expressed in the 
Ward identities. In particular, these Green functions can, 
in principle, be calculated for any value of squared mo- 
menta even though ChPT is set up only for a low-energy 
description. For the discussion of tttt scattering one con- 
siders the four-point function [?] 



Gabcd / \ 
PPPPy-'^ai Xfn Xc, Xd) 



{mp''{^a)---p\xd)m (5) 



with the pseudoscalar quark density — iqj^r'^q. In 
order so see that Eq. ^ can indeed be related to tttt scat- 
tering, we investigate the matrix element of evaluated 
between a single-pion state and the vacuum [?] : 



(0|F'^(0)|7r''(g)) =6''''G^ 



(6) 



The coupling of an external pseudoscalar source p to the 
Goldstone bosons is given by 



C,^t = i^TvipU^ - Up) 
f 2BFt^p'^tt°' , 



(7) 



where the first and second lines refer to the a and expo- 
nential parameterizations of U, respectively. From Eq. ^ 
we obtain = 2BF^ independent of the parameteriza- 
tion used which, since the pion is on-shell, is a consequence 
of the equivalence theorem [37] . As a consistency check, 
let us verify the PCAC relation from the QCD Lagrangian 

d^.A''"- = miq-i^T°-q = mP'' , 

evaluated between a single-pion state and the vacuum. 
The axial- vector current matrix element obtained from C2 
reads 

(0|A'^'"(a;)|7r^(q)) ^g^F^e-*«■^(5'^^ (8) 
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Taking the divergence implies M^F-^ = ttlGt^. In other 
words, 

, mP'^(x) 



M^F, 

can serve as a so-called interpolating pion field in the LSZ 
reduction formula. Using Eq. ([5]), the reduction formula 
relates the ^-matrix element of tttt scattering to the QCD 
Green function involving four pseudoscalar densities 

S,.= [0 {pl~Ml)--.{pl^Ml) 

X I d Xa ' • • d Xd e ^ • ■ ■ e ^*^ ppppiXa^ Xc^ ^d)- 



Using translational invariance, let us define the momen- 
tum space Green function as 

{2Tr)'^S'^{Pa +Pb + Pc + Pd)Fp^ppp{Pa,Pb,Pc,Pd) = 

j4_ —iPa-^a —iPh-^b —iPc-Xc —iPHXij 

d Xad Xbd Xcd Xiie ' e ' e e 

xGpppp^Xa, Xb, Xc, X(i), 

where we define all momenta as incoming. The usual rela- 
tion between the 5* matrix and the T matrix, S = I + iT , 
implies for the T-matrix element {f\T\i) = (27r)^(5'*(Py — 
Pi)Tfi and, finally, for M = iTfi: 



M = 



1 

Gi 



n 



k—a.b.c.d P^ 



lim (pi 



Ml) 



rpabcd 
^ PPPP 



(10) 



with FI^Ypp = Ff^l,p{pa,Pb, -Pc, -Pd)- We will now de- 
termine the Green function F^^ppp using the a and expo- 
nential parameterizations for U . In the first parameteriza- 
tion we only obtain a linear coupling between the external 
pseudoscalar field and the pion field [see Eq. ([7])] so that 
only one Feynman diagram contributes 



^PPPP — l^-D-f^Trj 



Pi - Pi - Mi 



(11) 



The Green function Fp''ppp depends on six independent 
Lorentz scalars which can be chosen as the squared invari- 
ant momenta p^ and the three Mandelstam variables s, t, 
and u satisfying the constraint s + i + u = ^j,p^. 

Using the second parameterization we will obtain a 
contribution which is of the same form as Eq. (fTTj) but with 
AA2^ replaced by M.\'^ . Clearly, this is not yet the same 
result as Eq. (fTT|) because of the terms proportional to Ak 
in However, in this parameterization the external 

pseudoscalar field also couples to three pion fields [see Eq. 
([7])], resulting in four additional contributions 



AaF^ 



bed 
PPPP 



A Tpabcd 
^d^ PPPP 



with 



AaF^ypp{pa,Pb,-Pc,-Pd) 

= i2BF^Y- 



Pl - Mi pI - Mi 



iA„ 



3Fi 



(12) 



and analogous expressions for the remaining AFs. In total, 
we find a complete cancelation with those terms propor- 
tional to Ak (in the second parameterization) and the end 
result is identical with Eq. Finally, using G-^ — IBF-j^ 
and inserting the result of Eq. ([TT]) into Eq. ((TU]) we obtain 
the same scattering amplitude as in the "direct" calcula- 
tion of A^l^ and A^2* evaluated for on-shell pions. 

This example serves as an illustration that the method 
of Gasser and Leutwyler generates unique results for the 
Green functions of QCD for arbitrary four-momenta. There 
is no ambiguity resulting from the choice of variables used 
to parameterize the matrix U in the effective Lagrangian. 
These Green functions can be evaluated for arbitrary (but 
small) four-momenta. Using the reduction formalism, on- 
shell matrix elements such as the tttt scattering amplitude 
can be calculated from the QCD Green functions. The re- 
sult for the TTTT scattering amplitude as derived from Eq. 
([TU)) agrees with the "direct" calculation of the on-shell 
matrix elements of Ai2^ and Ai'^'^- On the other hand, 
the Feynman rules of 7^2^ and A^2"^ when taken off shell, 
have to be considered as intermediate building blocks only 
and thus need not be unique. 



3.2 Model (in)dependence of pole positions and 
Breit-Wigner parameters 

A popular definition of masses of unstable particles cor- 
responding to a (relativistic) Breit-Wigner formula makes 
use of the zero of the real part of the inverse propagator. 
It has been shown that such a definition leads to field- 
redefinition and gauge-parameter dependence of the mass 
starting at two-loop order [38 , 3 9[|40[HIll42ll43ll44| . In con- 
trast, defining the mass and width in terms of the complex- 
valued position of the pole of the propagator leads to both 
field-redefinition and gauge-parameter independence. 

As the baryon resonances are thought to be described 
by QCD, with the progress of lattice techniques and, es- 
pecially, the low-energy effective theories (EFT) of QCD 
(see, e.g., [45 l l35 l l46 l [47 l l48 l l49 l [50] and references therein) 
the question of defining baryon resonance masses becomes 
important. Here we examine this issue for the A reso- 
nance. As discussed in Ref. [7], the conventional reso- 
nance mass and width determined from generalized Breit- 
Wigner formulas have problems regarding their relation 
to S-matrix theory and suffer from a strong model de- 
pendence. Here, we will show that these parameters, in 
addition, depend on the field-redefinition parameter in a 
low-energy EFT of QCD. 

For simplicity we ignore isospin and consider an EFT 
of a single nucleon, pion, and A resonance. Defining 

Af_,i, = -{i^~mA)g^,.+i {l^di,+-fi,dfj,)-ijf^^ji,~mAjf_,^^, 

the free Lagrangian is given by 



£0 = V A^, r jun)^ -f - a^TT^'^Tr . (13) 



Here, the vector-spinor describes the A in the Rarita- 
Schwinger formalism [51j . ^ stands for the nucleon field 
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with mass toat, and tt represents the pion field which we 
take massless to simplify the calculations. The interaction 
terms have the form 

Ant = 9 d-'TT i^^' {g^, - 7^7,) ir + H.c. + • • • , (14) 

where the ellipsis refers to an infinite number of interac- 
tion terms which are present in the EFT. These terms 
also include all counter-terms which take care of diver- 
gences appearing in our calculations. Although our results 
are renormalization scheme independent, for simplicity we 
use the dimensional regularization with the minimal sub- 
traction scheme. 

Let us consider the field transformation 

ipt' ^^pt' +^d^'Tr§ , ^iIj" +^d''Tr9 , (15) 

where ^ is an arbitrary real parameter. When inserted into 
the Lagrangians of Eqs. (fT3|) and (fT4| . the field redefinition 
generates additional interaction terms. The terms linear in 
^ are given by 

C^dd int = e § A^, r+i d-'n V If . (16) 

Because of the equivalence theorem physical quantities 
cannot depend on the field redefinition parameter ^. The 
complex-valued position of the pole of the A propagator 
does not depend on ^. In contrast, the mass and width 
defined via (the zero of) the real and imaginary parts of 
the inverse propagator depend on ^ at two-loop order. 

The dressed propagator of the Z\ in n space-time di- 
mensions can be written as 

_ J _ tjL _ P^7" - 7^P" _ (n-2)p^p- - 
n — 1 (n — l)m/i {n — 1)?ti^ 

X — ; h pole- free terms , (17) 

where we parameterize the self-energy of the A as 

Si{p^)g^'' + ^2(p')7^7'^ + ^3(p')p^7"^ + ^4(p')7V 
+S^{p^)py' + S^{p')Tl>g'''' + ^7(^2)^7^'^- 

+i:8(p2)^P^7" + ^9(p')#7>" + ^io(p')|^P>^ (18) 

The complex pole z of the A propagator is obtained by 
solving the equation 

z-mA- Si{z^)- zSf,{z^) = Q. (19) 

The pole mass is defined as the real part of z. 

On the other hand, the mass ma and width F of the A 
resonance are often determined from the real and imagi- 
nary parts of the inverse propagator (corresponding to the 
Breit-Wigner parametrization) , i.e., 

rri]^ — niA — Re Si{mj^) — rriji Re Z'g(TO^) = , 

r = -2ImZ'i(m^) ~2mRlmSe{mj^) ■ (20) 

We have calculated the A mass using both definitions and 
analyzed their ^ dependence to first order (for details see 
Ref. [52]). 




(c) (d) 

Fig. 5. A self-energy diagrams. Solid, dashed, and double 
lines correspond to nucleon, pion, and A, respectively. 

The A self-energy at one loop-order is given by the 
diagram in Fig. [5] (a) . The two-loop contributions to the A 
self-energy are given in Fig. [5](b) - (d). We are interested 
in terms linear in ^. 

To find the pole of the propagator we insert its loop 
expansion 

z = rriA + Siz + 62Z + • ■ • (21) 

in Eq. (|19p and solve the resulting equation order by order. 

The one-loop diagram results in the ^-independent ex- 
pression for 6iz. Calculating diagram (b) and (c) we find 
that they give vanishing contributions. The ^-dependent 
contributions in S2Z, generated by the one- loop diagram 
and by diagram (d) exactly cancel each other leading to 
the ^-independent pole of the propagator. 

We perform the same analysis inserting the loop ex- 
pansion of TTlR, 

mj^ — rriA + Sim + + • • • , (22) 

in Eq. (|20|) . For Smi the one- loop diagram gives a ^- 
independent expression. On the contrary, the ^-dependent 
contributions in 621^, generated by the one-loop diagram 
and by diagram (d) do not cancel, thus leading to a ^- 
dependent mass mn. An analogous result holds for the 
width r obtained from Eq. ((20)) . 

To conclude, we addressed the issue of defining the 
mass and width of the A resonance in the framework of a 
low-energy EFT of QCD. In general, the scattering ampli- 
tude of a resonant channel can be presented as a sum of the 
resonant contribution expressed in terms of the dressed 
propagator of the resonance and the background contribu- 
tion. The resonant contribution defines the Breit-Wigner 
parameters through the real and imaginary parts of the 
inverse (dressed) propagator. The resonant part and the 
background separately depend on the chosen field vari- 
ables, while the sum is of course independent of this choice. 
We have performed a particular field transformation with 
an arbitrary parameter ^ in the effective Lagrangian. In a 
two-loop calculation we have demonstrated that the mass 
and width of the A resonance determined from the real 
and imaginary parts of the inverse propagator depend on 
the choice of field variables. On the other hand, the com- 
plex pole of the full propagator does not depend on the 
field transformation parameter ^. 

The above conclusions are not restricted to the con- 
sidered toy model or EFT in general. Rather, our results 



S. Capstick et el: The physical meaning of scattering matrix singularities 



9 



are a manifestation of the general feature that the (rel- 
ativistic) Breit-Wigner parametrization leads to model- 
and process-dependent masses and widths of resonances. 
Although in some cases (like the A resonance) the back- 
ground has small numerical effect on the Breit-Wigner 
mass, still the pole mass and the width should be consid- 
ered preferable as these are free of conceptual ambiguities. 



F15 transverse helicity amplitudes 



4. 

4 



hypercentral CQM 
Photon point (PDG) 



A,,2Jlab05 
,2 tiypercentral CQM ■ 
Photon point (PDG) 
A3,2 Old data 

.u-ib 05 



4 Bare and dressed quantities within a well 
defined model 

4.1 Longitudinal and transverse helicity amplitudes of 
nucleon resonances in a constituent quark model - 
bare vs dressed resonance couplings 

Many models have been built and applied to the descrip- 
tion of hadron properties. An important role is played 
by Constituent Quark Models (CQM), in which quarks 
are considered as effective degrees of freedom. There are 
many versions of CQM, which differ according to the cho- 
sen quark dynamics: h.o. and three-body force [53l, alge- 
braic [54|, hypercentral [55], Goldstone Boson Exchange 
[55] . instanton [S7]- Here we report some results of the 
hypercentral CQM (hCQM) [55] on the longitudinal and 
transverse helicity amplitudes of the nucleon resonances. 
In this model, the quark interaction is assumed to be given 
by a hypercentral potential 



-50 - 



1 V- 



V(x) = — r/.T + 



= (23) 



where x is the hyperradius expressed in terms of the in- 
ternal Jacobi coordinates p and A. A coulomb-like plus 
linear confinement form of the potential is supported by 
recent lattice QCD evaluations of the quark-antiquark po- 
tential [58J and in this sense Eq. (|23)) can be considered as 
the hypercentral approximation of a two-body Cornell-like 
potential. The model interaction is completed by adding a 
standard spin dependent hyperfinc interaction H^yp [53 , 
in order to reproduce the splittings within the SU (6) mul- 
tiplets. The few free parameters (a, t and the strength of 
Hhyp) are fitted to the spectrum and the model is then 
applied to calculate (i.e. to predict) various properties of 
hadrons: the photocouplings |59], the transverse helicity 
amplitudes for negative parity resonances [60j , the elastic 
form factors [61], the longitudinal and transverse helicity 
amplitudes of all the main resonances [62j . 

It is interesting to analyze in a systematic way the 
behaviour of the helicity amplitudes in comparison 
with the existing data. Fig. [6| shows the results for the 
transverse helicity amplitudes of the F15 resonance, results 
which are typical for a J > 1/2 state [60''62' : the medium- 
high behaviour is quite well reproduced, showing that 
the hypercoulomb part of the interaction 1/x gives a fair 
account of the short range, while at low there is a 
lack of strength, particularly for the amplitude. For 
J = 1/2 states [59ll62] . there are some minor problems in 
the low region, but for the rest the agreement with 
data is satisfactory. Major problems are present for the 
Roper resonance [62l , a fact that may support the idea of 



1 2 3 4 6 

(GeV^) 

Fig. 6. The transverse helicity amplitudes for the F15 reso- 
nance obtained with the hCQM [62], compared with the ex- 
perimental data, taken by an old compilation [63j, recent JLab 
experiments |64] and PDG [65) . 



a particular status of the radial excitations of the nucleon. 
Discrepancies are present also for the Z\— resonance |62) . 
a feature which is typical of all CQMs; it is well known 
that the quark model, while reproducing quite well the 
baryon magnetic moments, fails in the case of the N — A 
transition magnetic moment. Taking into account the fact 
that the proton radius, calculated with the wave functions 
corresponding to the potential of Eq. ([23)) . turns out to be 
about 0.5 fm, the emerging picture is that of a small quark 
core surrounded by an external region, which is probably 
dominated by dynamical effects not present in the CQM, 
that is sea-quark or meson cloud effects ^60^ ■ 

These considerations are relevant in connection with 
the issue of bare vs. dressed quantities. One should not 
forget that the separation between bare and dressed quan- 
tities is meaningful within a definite theoretical approach. 
In CQM calculations the aim is not a fit but the descrip- 
tion of observables, which in principle are dressed quanti- 
ties (like baryon masses, magnetic moments, helicity am- 
plitudes, etc). In any case the identification of quark re- 
sults with bare quantities is questionable in view of the 
fact that CQs have a mass and some dressing is implicitly 
taken into account. However a consistent and systematic 
CQM approach may be helpful in order to put in evidence 
explicit dressing effects. 

These effects have been recently calculated by means of 
a dynamical model [66] . The meson cloud contribution to 
various helicity amplitudes has been calculated and com- 
pared with the hCQM predictions [67] . The two contri- 
butions cannot be added, since they are calculated within 
different frameworks, however it is interesting to note that 
the meson cloud contribution is relevant at low and in 
most cases it is important where the hCQM prediction 
underestimates the data. An example of this situation is 
given by the longitudinal and transverse helicity ampli- 
tudes for the Z\— excitation [67]. The case of the S'1/2 am- 
plitude is particularly interesting (see Fig. [7|) : the hCQM 
is almost vanishing and the meson cloud contribution ac- 
counts for practically the whole strength. 
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Fig. 7. The dependence of the N —* A longitudinal he- 
licity amplitude 51/2: superglobal fit performed with MAID 
[68] (solid curve), predictions of the hypercentral Constituent 
Quark Model [55ll67ll62] (dashed curve), pion cloud contribu- 
tions calculated with the Mainz dynamical model [66] (dotted 
curve). The data points at finite are the results of single-Q^ 
fits [67] on recent data quoted in ref. [67] . 



One should however be aware of some problems, both 
phenomenological and theoretical. From the phenomeno- 
logical point of view, there is the problem of the sign of 
the helicity amplitudes, which is actually extracted from 
the meson electroproduction amplitude, the fact that the 
PDG photon points are often non consistent and the need 
of new and systematic data. The main theoretical prob- 
lems are connected with the inclusion of relativity. The 
kinematic relativistic corrections seem to be not impor- 
tant for the helicity amplitudes f75', however relativity 
should be consistently included both in the (unquenched) 
CQM states and in the transition operators, leading to the 
possibility of quark pair terms in the electromagnetic cur- 
rent. In any case, the unquenching of the CQM is expected 
to produce a substantial improvement in the theoretical 
description of baryon properties. In particular, it will be 
possible to calculate simultaneously the electromagnetic 
processes and the strong decays and the baryons reso- 
nances will acquire a non zero width through the coupling 
to the continuum part of the spectrum. 



The problem is how to introduce dressing in the calcu- 
lations. One way is to adopt a hadronic approach: meson 
and baryons (nucleon and nucleon resonances) are the rel- 
evant degrees of freedom and the dynamics is given by 
meson-baryon interactions. This is certainly a consistent 
approach which has been used with success by various 
groups with different techniques (see e.g. [68,69,70])- An- 
other possibility is given by the so called hybrid models, 
where the baryons are considered as three-quark states 
surrounded by a pion cloud and a direct quark-meson cou- 
pling is introduced. In this way the electromagnetic exci- 
tation acquires contributions also from the meson cloud. 
This approach is very useful for preliminary calculations 
(see ref. 1711 ). however a more promising method is pro- 
vided by the inclusion of dressing mechanisms directly at 
the quark level. This means in particular the inclusion of 
higher Fock components in the baryon state: 



{I'd > = I'Sq \qqq > + ^3q qq I 99 > 



(24) 



and implies the necessity of unquenching the quark model, 
as discussed in section 2.1. For the case of mesons, there 
are pioneering works by Geiger and Isgur [72j . where the 
qq pair creation mechanism is introduced at the micro- 
scopic level within a string model. In the case of baryons, 
the problem is more complicated and has been recently 
treated performing the sum over intermediate quark loops 
by means of group theoretic methods [73]. This approach 
has been applied to the determination of the strange con- 
tent of the proton [71j with good results. In this way we 
have at our disposal a promising method for obtaining 
an unquenched, that is dressed, formulation of the CQM. 
The systematic calculation of baryon properties, such as 
transition amplitudes (but also elastic form factors and 
structure functions) in an unquenched CQM will supply 
a set of dressed quantities to be compared directly with 
data and will allow to understand where meson cloud or 
(better) qq effects are important. 



4.2 Nucleon Resonances and Hadron Structure 
Calculations 

4.2.1 What are the nucleon resonances ? 

To answer this question, it is useful to recall some text- 
book (For example, see Refs. fTSllTTlFfH] ) definition of reso- 
nances. Phenomenologically, a resonance (i?) is identified 
with a peak in a plot of the reaction cross section as func- 
tion of the collision energy E or invariant mass W . At en- 
ergies near the peak position, one can fit such data near 
the peak by 



R,a\ 



\E 



R,b\ 



iW-mR)^ + \^ 



2 



(25) 



where p{W) is an appropriate phase space factor, mn the 
position of the peak, and Fq the width of the peak. The 
expression Eq. ([25|) has the same function form of the de- 
cay rate of an unstable system with a mass ~ niR and a 
life time r/j ~ I/Iq. It is thus natural to interpret that the 
cross section Eq.(|25]) is due to the excitation of an unsta- 
ble system during the collision. How this unstable system 
is formed from the entrance and exit channels is a dynam- 
ical question which can only be answered by modeling the 
reaction mechanisms and the internal structure of all par- 
ticles involved. We will address this question within the 
Hamiltonian formulation of the problem. This is rather 
different from the S-matrix approach. 

Within the Hamiltonian formulation, there are two 
ways to derive the expression Eq. p5|) depending on the 
structure of the excited unstable system. Let us first con- 
sider the one defined in Feshbach's textbook (page 23 of 
Ref.JSj). It can be stated as the following : 

A resonance is formed in a process that the in- 
cident projectile completely lose its identity, amal- 
gamating with the target system to form a com- 
pound state. Namely, the evolution of the v^rhole 
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system can not be defined in terms of the motion 
of the projectile and its transmutation. 



The expression Eg. ([25)1 corresponding to this defini- 
tion of resonances can be formulated by assuming that 
the Hamiltonian of the system has the following form 



with 



(26) 



(27) 



where Fn^a defines the decay of an unstable system R 
with a mass Mq into channel a. The reaction amplitude is 
defined by 



T{E) ^ H' + H'- 



1 



H + ie 



-H' 



(28) 



^From Eqs. ([?7)l - (l^5|) . it is straightforward to see that the 
reaction cross section for b ^ a can be written as 



(29) 



with 



rR,a{ka)rR^b{h) 



W -Mo- S{W) 
where ka is the on-shell momentum of channel a and 

1 



(30) 



SRiw) = J2<R\rl^- 



'W-Ho + ie 
We can cast Eq. (|30|) into 



W-MR{W)+^^^^ 



rR^a\R> (31) 



(32) 



where 



Mr{W) = M" + Re{SR{W)) 

r*°*(w) 



(33) 
(34) 



An unstable system is formed during the colli- 
sion by an attractive force between the interacting 
particles which do not lose their identities. 

The simplest parameterization of an attractive force is 
a separable potential 



The solution of Eq. ((28|) then become 

9*ik)gik) 



where 



T{W) 



z{W) =< g\ 



C-z{W) 



1 



W -Hn + ie 



-\9> 



(37) 



(38) 



(39) 



If the parameters of H' are chosen such that C—Re{z(W)) - 
on the physical world W — > Wq where Wq is a real num- 
ber, we can expand 

C -z{W) = (40) 
[ C - R{Wo) - R'{Wo){W -Wo) + ---]- iI{W) 
- -R'iWo) 

where R{Wo) = Re{z{Wo)), /(VK) = Im{z(W)), and 
R'iWo) = dRe{z{W))/dW\w=Wo- We then can write at 
W -^Wo 



T {W Wq) 



(41) 



-9l(ko)R7W^9b{ko) 



W - [Wq + (RiWo) - C)] + i p/f^^i I{Wo) 



R'iWo) 



The above expression can give resonant cross section Eq.(|25l 
if the parameters of 9 and C can be chosen to satisfy 



mR = Wo 

ro 



1 



1 



2 R'{Wo 



R'iWo 

-nwo) 



■{R{Wo)-C) 



(42) 
(43) 



By using Eqs.([29l) and to fit the expression Eq.([25|. 
the parameters of the Hamiltonian are then related to the 
data by the following relations 

mR = Mr{W = tur) = M° + Re{SR{mR)) (35) 
To = r'^^iW = niR) = -2ImiIJR{mR)) (36) 



Eqs. (f35|) - ([36|) then allow us to use the experimental values 
niR and /q to extract the property of the unstable system, 
specified by Mo and FR^a of the Hamiltonian, through the 
evaluation of Eqs.([3T |l -([33 | at energies near W = itir. 

The second mechanism which can also yield a cross 
section of the form of Eq. ([25|l is : 



4.2.2 Dynamical Models for investigating Nucleon 
Resonances 

^From the above two examples, we see that the resonant 
cross section Ea. (l25p can correspond to two very differ- 
ent internal structure of the excited unstable system. The 
nucleon resonances we are interested in correspond to the 
unstable systems defined by the Hamiltonian Eq. (|?7)l . For 
the meson-nucleon reactions, such unstable systems are 
due to the excitation of the quark-gluon substructure of 
the nucleon. 

In reality, the situation is much more complicated. In 
the reactions involving composite systems, such as atoms. 
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nuclei and nucleons, the excitations of resonances always 
involve non-resonant direct interactions. For example, the 
non-resonant interactions in pion-nucleon scattering could 
be due to the exchange of p meson. The reaction formula- 
tion for analyzing such reactions is well presented in Fes- 
hbach's textbook[78]. We now briefly describe how such a 
formulation can be used to investigate nucleon resonances 
in meson-nucleon reactions. 

The starting point is to divide the Hilbert space into a 
P space for the entrance and exit channels and Q for the 
rest. One can cast the equation of motion in the P-space 
as 



{E - H,ff)P^ = 0, 



where 



^eff ~ Hpp + HpQ 



Ei+) - Hr 



(44) 



-Hqp . (45) 



Here E'^^^ = E + ie specifies the boundary condition and 
we have defined projected operator Hpp = PHP, Hpq = 
PHQ and Hqq = QHQ. Now consider the eigenstates of 
Hqq which can be discrete bound or unbound <l>e,a 
states 



(46) 
(47) 



with 



<^s\^s' > ^Ss^,, (48) 
< -Pie, a)\^ie', a') > = (5„,„.,5(e - e') (49) 

We then expand 
^eff ^Hpp = 

HpQ^s >< ^sHqp 



+ 



da / de 



E-e, 

HpQ^{e,a) >< ^{e,a)HQp 

EW - e 



(50) 



One can see from the above equation that rapid energy 
dependence of ife// will occur as the energy approaches 
one of the bound state energy e^. This is the origin of 
rapid energy-dependence of the cross sections. As shown in 
Feshbach's book (page 158-162), the amplitude at E' ~ 
can be written as 



< X^-MHpq\<Ps >< <P.\HoPY^+'^ > 



E - es- < <Ps\Wqq\<Ps > 



with 



where 



Wqq = H 



1 



QP 



E(+) - Hpp 



-H 



PQ 



(51) 



(52) 



Hpp = Hpp 



da I , Hp^^i^,a)^><^{e,a)H^p ^^^^ 



and x^^"* are the solutions of 



{E - Hpp)x^+^ = 
(E ~ Hpp)x''-> = 



(54) 
(55) 



In this formulation, one bound state of Hqq will corre- 
spond to one resonance. Namely, one can predict whether 
a resonance can appear in a particular partition of Hilbert 
space by examining whether bound states can be gener- 
ated from the Hamiltonian when the coupling with the 
states P are turned off. 

We now point out that the dynamical model developed 
in Ref.[8] (MSL model) for investigating meson-baryon 
reactions are completely consistent with the formulation 
given in Eqs. ([5T|l - ([551) . To see this, one just make the fol- 
lowing identifications: 

— P space contains reaction channels 

MB = ttN, jN, -qN, ttA, pN, aN and ttttN 

P = Y^ \MB >< MB\ + \t:ttN >< t:ttN\ (56) 

MB 

— Hqq describes the internal structure of the bare N* 
states 



Hqq\N* > ^ M^N* > 



(57) 



— Hpp defines the non-resonant meson-baryon interac- 
tions 



Hpp = J2 > yruB+p'^ + y/mB + k-^] < MB\ 

VmB,M'B' + ''^[vmB.ttttN + ^TnrN,MB] 
MB 

(58) 



MB 



MB,M'B' 



— Hqp defines the coupling of the internal structure of 
N* with the reaction channels 

Hqp = E[*-X] ^N',MB + Pn* ] (59) 
N* MB 

With some inspections, one can see that equations 
presented in the section 3 of Ref.[8] (MSL model) are 
completely equivalent to Eqs. ((5T|) -([55 l) . If we set MB = 
7tN,^N and Q = \A >< A\, and neglect ttttN channel, 
we then obtain the formulation of the SL model [11] . 



4.2.3 Relations with Hadron Structure Calculations 



We now note that Cs and <Ps in Ea. p6|) and (|5ip relate 
the structure calculations for the unstable systems in the 
Q-space and the reaction amplitudes in the P-space. In 
the MSL formulation, these are the bare mass M]^, and 

wave function |A^* > of the discrete bound states defined 
by Eq. ([57|) . These bare N* states can be considered as the 
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excited states of the nucleon if its coupling with the reac- 
tion channels is turned off. It is therefore natural to spec- 
ulate that the bare N* states in the SL and MSL models 
correspond to the predictions from a hadron model with 
confinement force, such as the well-developed constituent 
quark model with gluon-exchange interactions. This was 
first noticed in the SL model in 1996. The idea was later 
pursued in Ref. IJ^ in 2000 in an attempt to directly cal- 
culate ttN scattering amplitude in 6*11 up to = 2 GeV 
starting from several constituent quark models. In the 
later N{e,e'N) analysis 12,79 based on SL model, the 
bare jN A (1232) form factors were also found to 
be close to the constituent quark model predictions. To 
consider constituent quark models with meson-exchange 
residual interactions, the SL and MSL models must be 
modified to account for the contribution due to the contin- 
uum in the Q-space; namely the effects due to the second 
term in the right-hand-side of Eq. (f53|) . 

It is unlikely that the Lattice QCD calculation (LQCD) 
can account for the channel coupling effects and unitarity 
conditions, which are the essential elements of a dynamical 
coupled-channel analysis, rigorously in the near future. It 
is a challenging problem to relate the LQCD calculations 
to the information which can be extracted from the full 
solution Eq. (|5ip of a dynamical coupled-channel analysis. 

One possibility is to perform a LQCD calculation which 
defines a Hqq of a dynamical coupled-channel analysis. At 
the present time, perhaps the predictions from a quenched 
LQCD with heavy quark mass and no chiral extrapola- 
tion correspond to the bare parameters resulted from the 
dynamical coupled-channel analysis being performed at 
EBAC. This is based on the argument that the quark-loop 
contributions are suppressed at heavy quark limit and the 
LQCD mainly accounts for the gluonic interactions which 
are not in the P-space of MSL model. 



5 Conclusion 

This work was motivated by problems relating the reliable 
results from partial-wave and amplitude analysis, which 
are the parameters of dressed scattering matrix singulari- 
ties, and the results of quark models, which are usually 
given as the properties of bare resonances. Undressing 
dressed scattering matrix singularities in coupled-channel 
models involves model-dependent hadronic mass shifts, 
which arise from the unmeasurability of off-shell effects 
accompanying the dressing procedure. It is legitimate to 
extract bare quantities in coupled-channel models within 
the framework of a well defined model, but their inter- 
pretation requires keeping track of hadronic mass shifts 
produced by off-shell ambiguities. The best meeting point 
between quark model predictions and analyses of exper- 
imental data are dressed scattering matrix singularities, 
as although dressing (un- quenching) the quark model is 
complicated, it is in principle a solvable problem. This 
will require careful definition and checking of the proce- 
dures for extracting poles from energy-dependent partial 
waves or directly from partial-wave data. 



References 

1. S. Capstick and W. Roberts, Prog. Part. Nucl. Phys. 45, 
S241-S331 (2000), and references therein. 

2. H. W. Fearing, Phys. Rev. Lett. 81 (1998) 758. 

3. H. W. Fearing and S. Scherer, Phys. Rev. C 62 (2000) 
034003. 

4. S. Scherer and H. W. Fearing, Nucl. Phys. A 684 (2001) 
499. 

5. M. Gell-Mann and M. L. Goldberger, Phys. Rev. 96, 1433 
(1954); R. Haag, Phys. Rev. 112, 669 (1958); J. S. R. 
Chisholm, Nucl. Phys. 26, 469 (1961); S. Kamefuchi, L. 
ORaifeartaigh, and A. Salam, Nucl. Phys. 28, 529 (1961); 
S. Coleman, J. Wess, and B. Zumino, Phys. Rev. 177, 2239 
(1969); B. K. Jennings, Phys. Lett. 196, 307 (1987). 

6. G. Hohler, ttN Newslett. 10, 320 (1997); 14, 168 (1998), 
and references therein. 

7. G. Hohler, Against Breit-Wigner parameters — a pole-emic, 
in C. Caso et al. [Particle Data Group], Eur. Phys. J. C 3, 
624 (1998). 

8. A. Matsuyama, T. Sato and T. -S. H. Lee, Physics Report 
493, 193 (2007), and references therein. 

9. Guan Yen Chen, S. S. Kamalov, Shin Nan Yang, D. Drech- 
sel, and L. Tiator L., Phys. Rev. C 76, 035206 (2007); Tia- 
tor, D. Drechsel, S. Kamalov, M.M. Giannini, E. Santopinto, 
A. Vassallo, Eur.Phys.J. A19 (2004) 55-60, and references 
therein. 

10. I.R. Afnan, Acta Physica Polonica 8, 2397 (1998); A.W. 
Thomas, Adv. Nucl. Phys. 13, 1 (1984); Dae Sung Hwang 
and Do-Won Kim, Phys. Lett. B 601, 137 (2004); Yu. S. 
Kalashnikova, Phys. Rev. D 72, 034010 (2005); P.,]. Fink, 
Jr., G. H, R. H. Landau and J.W. Schnik, Phys. Rev. C 41 
(1990) 2720. 

11. T. Sato and T.-S. H. Lee, Phys.Rev., C54, 2660 (1996) 

12. T. Sato and T.-S. H. Lee Phys.Rev., C63,055201 (2001) 

13. T. Yoshimoto, T. Sato, T.-S. H. Lee, Phys. Rev. C61, 
065203 (2000) 

14. P. Geiger and N. Isgur, Phys. Rev. D 41 (1990) 1595; Phys. 
Rev. Lett. 67 (1991) 1066; Phys. Rev. D 44 (1991) 799; Phys. 
Rev. D 47 (1993) 5050. 

15. M. A. Pichowsky, S. Walawalkar and S. Capstick, Phys. 
Rev. D 60 (1999) 054030. 

16. E.S. Swanson, private communication. 

17. W. Blask, M. G. Huber and B. Metsch, Z. Phys. A 326 
(1987) 413. 

18. P. Zenczykowski, Ann. Phys. (NY) 169 (1986) 453. 

19. B. Silvestre- Brae and C. Gignoux, Phys. Rev. D 43 (1991) 
3699. 

20. Y. Fujiwara, Prog. Theor. Phys. 88 (1992) 933; Prog. 
Theor. Phys. 89 (1993) 455. 

21. N. A. Tornqvist and P. Zenczykowski, Phys. Rev. D 29 
(1984) 2139; Z. Phys. C 30 (1986) 83. 



22. D. Morel and S. Capstick, ar Xiv:nucl-th/0204014| PIN 
Newslett. 16 (2002) 232. 

23. S. Capstick and N. Isgur, Phys. Rev. D 34 (1986) 2809. 

24. S. Capstick and W. Roberts, Phys. Rev. D 49, 4570 (1994); 
Phys. Rev. D 47, 1994 (1993). 

25. Ted Barnes and E.S. Swanson, talk presented at the ECT* 
workshop "Highly Excited Hadrons" (July 2005, ECT*, 
Trento, Italy). 

26. Work done in collaboration with S. Krewald, U.-G. 
Meifiner, and A. Wirzba. 



14 



S. Capstick et el: The physical meaning of scattering matrix singularities 



27. R. Haag, Phys. Rev. 112 (1958) 669. H.J. Borchers, Nuovo 
Cimento 25 (1960) 270. D. Ruelle, Helv. Phys. Acta 35 
(1962) 34. S. Coleman, J. Wess and B. Zumino, Phys. Rev. 
177 (1969) 2239. C. G. Callan, S. Coleman, J. Wess and 
B. Zumino, Phys. Rev. 177 (1969) 2247. 

28. V. Lensky, V. Barn, J. Haidenbauer, C. Hanhart, 
A. E. Kudryavtsev and U. G. Meissner, Eur. Phys. J. A 27, 
37 (2006) arXiv:nucl-th/0511054 . 

29. E. Epelbaum, Prog. Part. Nucl. Phys. 57, 654 (2006) 

30. S. Weinberg, Phys. Rev. 130 (1963) 776; V. Barn, 
J. Haidenbauer, C. Hanhart, Yu. Kalashnikova and 
A. E. Kudryavtsev, Phys. Lett. B 586, 53 (2004). 

31. O. Krehl, C. Hanhart, S. Krewald and J. Speth, Phys. Rev. 
C 62, 025207 (2000) 

32. B. Julia-Diaz, T. S. Lee, A. Matsuyama and T. Sato, 
larXiv:0704.16T5l [nucl-th] . 

33. S. Ceci, A. Svarc and B. Zauncr, Phys. Rev. Lett. 97 (2006) 
062002 

34. S. Scherer and H. W. Fearing, Phys. Rev. C 51 (1995) 359. 

35. J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

36. S. Scherer, Adv. Nucl. Phys. 27 (2003) 277. 

37. J. S. R. Chisholm, Nucl. Phys. 26 (1961) 469; S. Kame- 
fuchi, L. O'Raifeartaigh and A. Salam, Nucl. Phys. 28 (1961) 
529. 

38. S. Willenbrock and G. Valencia, Phys. Lett. B 247, 341 
(1990). 

39. G. Valencia and S. Willenbrock, Phys. Rev. D 42, 853 

(1990) . 

40. A. Sirlin, Phys. Rev. Lett. 67, 2127 (1991). 

41. A. Sirlin, Phys. Lett. B 267, 240 (1991). 

42. S. Willenbrock and G. Valencia, Phys. Lett. B 259, 373 

(1991) . 

43. J. Gegelia, G. Japaridze, A. Tkabladze, A. Khelashvili and 
K. Turashvili, in Quarks '92, Proceedings of the Interna- 
tional Seminar on Quarks, Zvenigorod, Russia, 1992, edited 
by D. Yu. Grigoriev, V. A. Matveev, V. A. Rubakov, and 
P. G. Tinyakov (World Scientific, River Edge, N. J., 1993), 
|ArXiv:hep-ph/9910527, 

44. P. Gambino and P. A. Grassi, Phys. Rev. D 62, 076002 
(2000). 

45. S. Weinberg, Physica A96, 327 (1979). 

46. J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B307, 
779 (1988). 

47. S. Scherer, Adv. Nucl. Phys. 27, 277 (2003); S. Scherer 
and M. R. Schindler, arXiv:hep-ph/0505265 

48. T. R. Hemmert, B. R. Holstein, and J. Kambor, J. Phys. 
G 24, 1831 (1998). 

49. C. Hacker, N. Wies, J. Gegelia, and S. Scherer, Phys. Rev. 
C 72, 055203 (2005). 

50. V. Pascalutsa, M. Vanderhaeghen, and S. N. Yang, Phys. 
Rept. 437, 125 (2007). 

51. W. Rarita and J. S. Schwinger, Phys. Rev. 60, 61 (1941). 

52. D. Djukanovic, J. Gegelia, and S. Scherer, Phys. Rev. D 
76, 037501 (2007). 

53. N. Isgur and G. Karl, Phys. Rev. D18, (1978) 4187; Phys. 
Rev. D 19, (1979) 2653; S. Capstick and N. Isgur, Phys. 
Rev. D 34, (1986) 2809; F. Cardarelh, E. Pace, G. Salme 
and S. Simula, Phys. Lett. B 357, (1995) 267. 

54. R. Bijker, F. lachello and A. Leviatan, Ann. Phys. (N.Y.) 
236, (1994) 69. 

55. M. Ferraris, M.M. Giannini, M. Pizzo, E. Santopinto and 
L. Tiator, Phys. Lett. B364, (1995) 231. 



56. L. Ya. Glozman and D.O. Riska, Phys. Rep. C268, (1996) 
263. 

57. U. L"oring, K. Kretzschmar, B. Ch. Metsch, H. R. Retry, 
Eur. Phys. J. AlO, (2001) 309; U. Loring, B. Ch. Metsch, H. 
R. Retry, Eur. Phys. J. AlO, (2001) 395; U. Loring, B.Ch. 
Metsch, and H. R. Retry, Eur. Phys. J. AlO, (2001) 447. 

58. G. Bah et al, Phys. Rev. D51, (1995) 5165; G. Bah, Phys. 
Rept. 343, (2001) 1. 

59. M. Aiello, M. Ferraris, M.M. Giannini, M. Pizzo and E. 
Santopinto, Phys. Lett. B387, (1996) 215. 

60. M. Aiello, M. M. Giannini, E. Santopinto, J. Phys. G: Nucl. 
Part. Phys. 24, (1998) 753. 

61. M. De Sanctis, E. Santopinto, M.M. Giarmini, Eur. Phys. 
J. Al, (1998) 187; M. De Sanctis, M.M. Giannini, L. Repetto, 
E. Santopinto, Phys. Rev. C62, (2000) 025208; M. De Sanc- 
tis, M.M. Giannini, E. Santopinto, A. Vassallo, Nucl. Phys. 
A 755, (2005) 294; Nucl.Phys. A 782, (2007) 57. 

62. M.M. Giannini, E. Santopinto, to be published. 

63. V. Burkert, private communication. 

64. I.G.Aznauryan et. al., Phys. Rev. C 72, (2005) 045201. 

65. W.-M. Yao et al. (PDG), J. Phys. G 33, (2006) 1. 

66. S. Kamalov, S.N. Yang, D. Drechsel, O. Hanstein, L. Tia- 
tor, Phys. Rev. C 64, (2001) 032201. 

67. L. Tiator, D. Drechsel, S. Kamalov, M.M. Giannini, E. 
Santopinto, A. Vassallo, Eur. Phys. J. A19, (2004) 55. 

68. D. Drechsel, O. Hanstein, S.S. Kamalov, L. Tiator, Nucl. 
Phys. A645, 145 (1999). 

69. S. Kamalov, S. N. Yang, D. Drechsel, O. Hanstein, L. Tia- 
tor, Phys. Rev. bf C 64, (2001) 032201. 

70. A. Matsuyama, T. Sato, T.-S. H. Lee, Phys. Rep. 439, 
(2007). 

71. D.Y. Chen, Y.B. Dong, M.M. Giannini, E. Santopinto, 
Nucl. Phys. A 782, (2007) 62. 

72. P. Geiger, N. Isgur, Phys. Rev. Lett. 67, (1991) 1066; Phys. 
Rev. D 44, (1991) 799; Phys. Rev. D 47, (1993) 5050; Phys. 
Rev. D 55, (1997) 299. 

73. E. Santopinto, R. Bijker, Proceedings of 11th Conference 
on Problems in Theoretical Nuclear Physics, Cortona, Italy, 
11-14 Oct 2006 (World Scientific, Singapore 2007) 331 (e- 
Print: hep-ph/0701227); invited Talk at the NSTAR2007 
Workshop and to be published. 

74. R. Bijker, E. Santopinto, Invited talk at 30th Sympo- 
sium on Nuclear Physics, Cocoyoc, Morelos, Mexico, 3-6 
Jul 2007 (e- Print: nucl-th/0703053) ; NSTAR2007 Workshop, 
contributed paper. 

75. M. De Sanctis, E. Santopinto, M.M. Giarmini, Eur. Phys. 
J. A2, (1998) 403. 

76. M. Goldberger and K. Watson, Chapter 8 of Collision 
Theory (Robert E. Krieger Publishing Company, Inc., 1975) 

77. A. Bohr and M. Mottelson, Appendix 3F of Nuclear 
Structure (W.A. Benjamin, Inc., 1969) 

78. Herman Feshbach, Chapters 1.4 and HI. 2 of Theoretical 
Nuclear Physics: Nuclear Reacttons (John Wiley & Sons, 
Inc., 1992) 

79. B. Julia-Diaz, T.-S. H. Lee, T. Sato, L.C. Smith Phys.Rev. 
C75, 015205 (2007) 



